Baryogenesis via Leptogenesis in an inhomogeneous universe 
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I. INTRODUCTION. 

Recent oscillation experiments [3, 0] have confirmed 
that the neutrinos have small but nonvanishing masses. 
An observation of neutrinolcss double beta decay by the 
future experiments 0] would hint toward the existence 
of very heavy right-handed Majorana neutrinos which 
generate naturally small neutrino masses via the see-saw 
mechanism. The existence of the heavy neutrinos can 
also explain [J] the observed baryon asymmetry of the 
universe since the three Sakharov conditions Q are eas- 
ily fulfilled in the standard model supplemented by the 
right-handed Majorana neutrinos. 

The see-saw and baryogenesis via leptogenesis mecha- 
nisms are very attractive from a theoretical point of view. 
Unfortunately, there is still no compelling evidence for 
the existence of the heavy right-handed Majorana neu- 
trinos. Apart from the neutrinolcss double beta decay 
experiments one can search for its signatures in accel- 
erator experiments or in lepton flavor violating decays 
@, 0] • However the corresponding amplitudes are highly 
suppressed by the large masses of the right-handed neu- 
trinos. 

During leptogenesis, processes with right-handed neu- 
trinos in initial, intermediate or final states are not sup- 
pressed. In a completely homogeneous and isotropic 
universe the outcome of baryogenesis via leptogenesis is 
a single number — the baryon-to-photon ratio, which 
is constant over space. In an inhomogeneous universe 
the primordial perturbations of the energy density and 
space-time metric modulate leptogenesis and the gen- 
erated baryon-to-photon ratio depends on the spatial 
coordinates. Thus instead of a single number we get a 
function of space coordinates which can, in principle, pro- 
vide us with some additional information on the model 
parameters. 
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A general perturbation with arbitrary fractional over- 
density in each component can be represented as a sum 
of adiabatic and isocurvature perturbations. For a purely 
isocurvature perturbation the sum of the fractional over- 
densities is zero. For a purely adiabatic perturbation 
the fractional overdensity in each matter component 
(baryons, photons, etc.) is the same and hence the 
baryon-to-photon ratio is constant over space. Conse- 
quently, a deviation of the baryon-to-photon ratio from 
its average value amounts to a baryon isocurvature per- 
turbation. We estimate here amplitude and scale depen- 
dence of the generated perturbations. 

Baryon isocurvature perturbations are known to affect 
primordial nucleosynthesis. The abundances of the pro- 
duced elements sensitively depend on the local baryon- 
to-photon ratio and the final average abundances can 
be considerably different from those in the homogeneous 
universe [1, 0, [T3, EH- They also affect the spectrum 
of the cosmic microwave background (CMB) radiation. 
Purely isocurvature perturbations produce a series of 
peaks in the CMB spectrum whose angular scales are 
approximately in the ratio 1 : 3 : 5 . . ., whereas purely 
adiabatic density perturbations produce peaks whose an- 
gular scales are in the ratio 1:2:3... [l2| ■ Current 
large-scale observations are consistent with the primor- 
dial perturbations being adiabatic. It is however possible 
that the upcoming experiments will find a small admix- 
ture of the isocurvature perturbations at scales smaller 
than 10 Mpc [HI, [3]; in particular the Planck satellite 
will be able to limit the amplitudes of isocurvature modes 
to less than 10% of the adiabatic mode [l||, HH. The 
isocurvature component modifies the sound velocit y in 
the plasma and the initial conditions for the CMB \u\ . 
and is expected to "smear out" the spectrum, i.e. change 
height and width of the peaks. 

The outline of the paper is as follows. In section HT1 we 
derive differential and integral forms of the Boltzmann 
equation and prove the gauge invariance of the latter 
one. The derived integral Boltzmann equation is used 
to analyze leptogenesis in section IIII1 where the gener- 
ated lepton asymmetry is related to the adiabatic energy 
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density and metric perturbations and possible velocity 
perturbation in the heavy neutrino component. Finally 
we summarize our results in section ITVl 



II. THE BOLTZMANN EQUATION 

The general form of the Boltzmann equation is [lE 
df 



dX 



C[f] 



(1) 



where A is the affine parameter along a geodesic and C 
the collision operator acting on the one-particle distri- 
bution function /. 

Since the Boltzmann equation treats particles as on- 
shell states, / is a function of space-time coordinates x a 
and three independent components of the particle mo- 
mentum Pi, or, alternatively, the particle kinetic energy 
E and unit momentum vector pi. Thus we can rewrite 
the left-hand side of fTJ) using the chain rule 



df df P i df 1 dE df 1 dp, _ C[f] 
dx° dx l P° dE P° dX dp, P° dX P° ' 



(2) 



The derivatives of E and pi can be calculated using the 
geodesic equation; the result depends on the space-time 
metric. 

The space-time metric of the early universe, which was 
almost homogeneous and isotropic, can be split into a 
background part 



g° = a 2 (77)diag(l, -1,-1,-1), 



(3) 



where rj is the conformal time, and a small perturbation 
which reads [TtJ 
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B-i 2 (ip6ij + Xij) 



(4) 



where a is the cosmic scale factor, 0, %j) and B arc scalar 
functions, and Xij defined by 



Xij 



did, - -^-OA ) £ 



(5) 



is traceless. Note that the definition of ip used here dif- 
fers by the ^dkdk£ term from that used in [l7|. This 
additional term will appear in the equation for the gauge 
transformation (|13aj) of ip later on. 

A nonzero shift function B means that comoving world- 
lines and worldlincs orthogonal to hypersurfaces of con- 
stant time are not collinear and we are dealing with a lo- 
cally nonorthogonal coordinate system [l^l • We avoid the 
associated unnecessary complications by setting B = 
which is fulfilled, for instance, in all synchronous gauges 
and the longitudinal gauge. 

Isotropy of the background space-time implies that the 
zero-order phase-space distribution functions are inde- 
pendent of the direction of the momentum vector, so 



that Jj- vanishes in the homogeneous universe defined 
by g® and is of first order in the inhomogeneous uni- 
verse defined by (jlj. Isotropy also implies that in the 
homogeneous universe the unit momentum vector is con- 
served, so that in the inhomogeneous universe -f± is of 
first order as well. Consequently, the product of the two 
is of second order and the last term on the left-hand side 
of @ can be neglected. Then the Boltzmann equation 
takes the form 

-,2 r 



df df P l _ df_ 
■C[f] 



P_ 

E 



H-ip-XiiPl 



+ <P,i PPi 




(6) 

The expression in the round brackets is (up to a sign) 
ptr4f calculated to linear order in the small perturba- 
tions of the metric, and H = ^ is the Hubble parameter. 
In the longitudinal gauge B = Xij = and the Boltzmann 
equation |6]) reverts to that derived in fl8l ]. 

As we are only interested in the total generated asym- 
metry, we integrate the Boltzmann equation ([6]) over the 
phase space. Integration of the first term on the left-hand 
side gives the derivative of the particle number density n 
with respect to the time coordinate x°. 

Since for a massive particle -pu = jjtt = V 1 , the second 
term on the left-hand side can be rewritten in the form 



df P 



d 



dx % 



d{nv l ) 
dx % 



(7) 



where v l denotes the macroscopic three-velocity of the 
gas. In the early universe the macroscopic three-velocity 
is a first-order quantity and therefore (up to second order 
corrections) is proportional to the corresponding macro- 



scopic four-velocity: v % ~ au l (recall that u 



Since Xij is °f first order, one can integrate the fifth 
term in ^ over the isotropic zero-order phase-space dis- 
tribution function. In this approximation the integration 
gives the same result for any i and therefore (recall that 
Xij is traceless) this term vanishes to first order. 

Analogously, since cf> is of first order, one can again in- 
tegrate the last term on the left-hand side of ([6]) over the 
phase-space using the zero-order phase-space distribu- 
tion function; the integral vanishes due to the isotropy of 
the distribution function. 

Upon integration by parts, the remaining terms take 
the form 



3(H — i')n = 



^=a 3 (l-3V0. (8) 



v -,g 3 dx° 

Using the fact that Xij is traceless it is straightforward 
to check that is the determinant of the spatial part of 
the metric. 

Collecting all the terms we obtain the integral form of 
the Boltzmann equation 



1 d(ny/-g 3 ) d(nu l ) 



dx° 



dx l 



E p 



(9) 
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Equation ([9]) is a partial differential equation which is 
rather difficult to solve. However, if we neglect some fur- 
ther second-order terms, it can be substantially simpli- 
fied. Since it* as well as the gradient of n are of first order, 
their product can be neglected, so that (nu l ).i ~ nu\. 
Next we introduce 



T 



a = V~93 ex P (/ u ]i adx° 



(10) 



As can be checked by direct substitution, equation (J9j) is 
equivalent (up to the second order terms) to the following 
equation for T 



1 dT 



cm 

E 



(11) 



where all the terms are evaluated at the same point in 
space. 

Equation pip is an ordinary differential equation with 
respect to the time coordinate. Any dependence of the 
generated lepton asymmetry on the space coordinates is 
"hidden" in the dependence of the metric and macro- 
scopic parameters on the space coordinates. This sub- 
stantial simplification is possible because we limit our 
analysis to linear perturbations here. 

From (jlip it is clear that T is conserved if the collision 
terms vanish; that is, in the inhomogeneous universe T 
replaces the particle number density per comoving vol- 
ume Y . 

Let us also note that equations © and (|11[) describe 
the development of the particle number density of a par- 
ticular species and can be used in a universe with an 
arbitrary mixture of adiabatic and isocurvature pertur- 
bations. 

The solutions of the Einstein equations give macro- 
scopic parameters like the temperature as functions of 
the time coordinate. Since the collision terms depend ex- 
plicitly on the inverse temperature z = 4lr (where Mi is 
mass of the lightest Majorana neutrino) j2CJ rather than 
on the time coordinate x°, it is convenient to rewrite (jlip 
as follows: 



C[f] 
E 



(12) 



The dot denotes differentiation with respect to the time 
coordinate, whereas the prime denotes differentiation 
with respect to z. 

The analysis of processes in the early universe is often 
complicated by the gauge choice. It is therefore impor- 
tant to check if the derived integral form of the Boltz- 
mann equation is invariant with respect to the gauge 
transformations which preserve B = 0. Let us consider 
a transformation from synchronous to longitudinal co- 
ordinates. Small perturbations of the metric and fluid 
velocity in the longitudinal gauge are related to those in 
the synchronous gauge by [17| 



~-d k dk£( s ), 



+ a di£. 



(13a) 
(13b) 



The same value of the proper time t (which is gauge 
invariant by definition) corresponds in the two gauges 
to two different values of the time coordinate: x® s ^ in 
the synchronous gauge and in the longitudinal one. 
Using the transformations laws for B and £ 

B (0 =B W +C + a#), £ { i)=£ {s) +C (14) 

(see [13 for more details) and taking into account that 
B is zero in both gauges and £ is zero in the longitudinal 
gauge, we obtain 



Ax° — x° s ) — £( s ). 



(15) 



Using equations ©, (|13ap and (fT5|) we find the following 
transformation law for the determinant of the spatial part 
of the metric 



-93 W (x°i)) = v / ^53( s) (a;( , s )) (1 + d k d k £ {s) ) . (16) 



Taking the divergence of (|13bp and multiplying by the 
scale factor we obtain 



(x° (l) ) = au l As) (x° {s) ) - d k d k £ 



(17) 



After integration over the time coordinate in (|10p and 
asymptotic expansion of the exponent, the second term 
on the right-hand side of (fl7|) cancels the analogous term 
in (|16p . so that a^{xn\) = <7( s ) (x^a). In other words, 
a is invariant under the gauge transformation. Three- 
dimensional scalars like the particle number density are 
gauge invariant as well: n^(x^) = "■( s )( a; ( ) s )) [HI- To- 
gether with the gauge invariance of a this implies the 
gauge invariance of T. 

Another important gauge invariant three-dimensional 
scalar is the energy density p 

P = T a = J.fV^gP°dP 1 dP 2 dP 3 = JfEdn p . (18) 

As can be shown by a direct calculation, for a gas of mass- 
less particles distributed according to the Boltzmann dis- 
tribution the leading contribution to the energy density 
due to a small macroscopic velocity u is proportional to 
u 2 , i.e. is of second order and can be neglected. Thus 
p = const. ■ X" 4 , and the temperature T (as well as the 
inverse temperature z) is gauge-invariant. 

Therefore the left-hand side of fi~2]) is gauge invariant 
since the derivative of one invariant quantity with respect 
to another invariant quantity is also invariant. 

On the right-hand side of (TT2|) the integral of the col- 
lision terms over the phase space is also invariant under 
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the gauge transformation 1 . Finally, since 



' dz 



(19) 



and the (inverse) temperature z as well as the proper 
time t are gauge invariant, the ratio (|19[) is also invari- 
ant. This completes our proof of the invariance of the 
integral Boltzmann equation under the gauge transfor- 
mation p3p . 

Instead of the transformation (fT5|) between syn- 
chronous and longitudinal gauge, we could have consid- 
ered a transformation between two different gauges which 
fulfill B = 0. The calculation would differ in only £ being 
now different from zero in both gauges and appearing also 
on the left-hand side of equations (fT6|) and (|T7|) . which 
however does not change the conclusion concerning the 
invariance of a. Therefore, we may conclude that equa- 
tion (|12p is invariant not only under the transformation 
between longitudinal and synchronous gauge, but under 
all gauge transformations that preserve the gauge condi- 
tion B = 0. 



III. INHOMOGENEOUS LEPTOGENESIS 

To calculate the lepton asymmetry generated in the 
decays of the heavy Majorana neutrinos one has to solve 
a coupled system of Boltzmann equations for the heavy 
neutrino and lepton number densities. 

To be able to perform this calculation analytically we 
neglect the contribution of the heavy neutrinos to the 
energy density of the universe. This is justified by the 
smallness of the associated number of effective massless 
degrees of freedom as compared to that of all the SM 
species. The two-body scattering processes and the fla- 
vor effects are also neglected for simplicity. In the sim- 
plest scenario of inflation isocurvature perturbations in 
the gas of the SM species are suppressed and we also ne- 
glect them. In particular, this implies that a is the same 
for leptons and photons so that T)l = ni/n 7 = T^/T 7 . 
The possible deviation of the macroscopic velocity of the 
heavy neutrino gas from that of leptons (i.e. velocity 
perturbation in the right-handed neutrino component) 
is parameterized by r u = ^ — 1 w J ' {u\ — u^) ;2 : adx°. 



A. Approximate analytical solution of the system 
of Boltzmann equations 

The collision terms are determined by the coupling of 
the Majorana neutrino to leptons and the Higgs 



C = N{f\^h) - \NMN C + h.c. 



(20) 



where I denotes the leptons, h = i^h^ is the charge 
conjugate Higgs doublet, and N are the components of 
physical Majorana neutrino field ^ = (N + N c ). 

Under the usual assumption that the phase-space dis- 
tribution function of the Majorana neutrino is propor- 
tional to the equilibrium one, the contribution of the de- 
cay (^ — > lb)) and the inverse d ecay [l\v — > \&) processes 
to the lepton asymmetry reads [201 



— ^-d£l p — (r^j 



,eq 



/ eq\ '"ilr 

e(ny - n^ 1 ) zqn e 



(21) 



where e is the usual CP asymmetry parameter and (Tx& 1 ) 
is the thermally averaged heavy neutrino decay width. 
Integrating over the equilibrium distribution function 
and neglecting a second-order contribution due to the 
small macroscopic velocity 2 we find 



= r, 



Ki{z) 



*i = 



(AAt) 



li 



8- 



(22) 



where K\ and K 2 are the modified Bcsscl functions. The 
equilibrium number densities are obtained by integra- 
tion of the corresponding Maxwell-Boltzmann distribu- 
tion functions over the phase space 



2 K 2 {z), 



cq 



AN- 



(23) 



where N = 3 is the number of generations. One should 
keep in mind that T refers to leptons and the Higgs; the 
right-handed neutrino is out of equilibrium and does not 
have a definite temperature. 

The same processes also contribute to the right-hand 
side of the Boltzmann equation for the Majorana number 
density [H 



E 



■dflp = -(r* x ) (n* - n^ q ) 



(24) 



In the homogeneous and isotropic universe dominated 
by radiation (see [l7| for more details) 



1 Perturbations of the space time metric may in principle lead to 
a change of decay widths and cross sections of the scattering 
processes. These effects, however, are not considered here. Fol- 
lowing common practice we assume that the widths and cross 
sections are given by the same expressions as those in vacuum. 
Then the perturbations only affect the thermally averaged decay 
widths and scattering cross sections through the modification of 
the temperature and the related modification of the equilibrium 
phase— space distribution functions. 



2 In the gas rest frame the thermally averaged decay width depends 
only on temperature. After a Lorentz boost to the initial refer- 
ence frame the decay width acquires the Lorentz factor, which 
differs from unity only by a second order term. 
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where (zo) Tq is the (inverse) background temperature 
and H = H(Mi) is value of the Hubble parameter at a 
temperature equal to the mass of the lightest Majorana 
neutrino. In the inhomogeneous universe we tentatively 
write 



f7 0cT- 3 (l + rv), 



/Soo 



£(l + r») (26) 



The formal solution of this equation satisfying the initial 
conditions A(0) = T^(0) = reads 



T L = -e (l+r u )T^ cq cxp (1 + r z ) lL dz") dz> . 

Jo 

(32) 

For z> 1 one can use the large-argument asymptotics of 
the Bessel functions K\ and K 2 , so that approximately 



where r a and r z are first-order gauge-invariant 3 func- 
tions. In what follows we will use r aqi and r a for the 
heavy neutrino and SM species, respectively. Introduc- 
ing finally k = Tqi 1 /'H, A = — T^ q and 



1d{z) = 



_ K^z) 



K 2 (zY 



lL (z) ee W K 2 (z) (27) 



we obtain the following system of Boltzmann equations 
for the lepton asymmetry and Majorana neutrino number 
density 



A' 



-kz (1 + r z )~f D A 



'y- cq 

1 * 



T' L = kz(1 + n) [e(l + r u ) lD A - 7L T L ] 



(28a) 
(28b) 



2tt 3 



exp(— z). (33) 



-r;\z)^[l + r^{z)-r^{z)] 



The arguments of the exponents in (|3"2"|) and (j3"3")l combine 
to a function which has a sharp peak at z = zj, where 
Zf is the solution of 



KZfy L (zf)[l + ri(z f )] = 1. 



(34) 



Consequently one can apply the method of steepest de- 
scent to estimate the asymptotic value of this integral. 
Approximating 7^ by its large-argument asymptotics 
and keeping only the terms of highest order in Zf we 
obtain 



where the relation r u = r^ — ru^ has been used. If v a , r a<Sl 
and r z vanish, the resulting system coincides with that in 
the homogeneous universe and the generated asymmetry 
depends only on e and k. 

In the weak washout regime (k <C 1) the deviation from 
equilibrium is large and washout processes play almost 



no role, so that T', 



-8(1 



Together with 



the initial conditions A(0) = T^(0) = this gives upon 
integration by parts and use of the zero-order solution of 
(I2~5al) . T%(z) = T*(0)exp(-Kz 2 /2), (see H) 



1 L 



l + r ff (0) + / r' u cxp {-nz 2 /2)dz 



(29) 



where the unperturbed value of the generated asymmetry 
in the weak washout regime is given by 



1 L 



2s 
^2" 



(30) 



In the strong washout regime (n 3> 1) the deviation 
from equilibrium is small, so that the final lepton asym- 
metry is independent of the initial conditions and deter- 
mined by the dynamics of leptogenesis. In this regime 
the density of Majorana neutrinos closely tracks its equi- 
librium value, which implies A' ~ (22| . Substitution of 
(f2"8"aj) into (|28b|) then gives 



T', 



-6(1 



kz (1 + r z )^ L Ti 



(31) 



-(Kz lL [l + r z ])' z ._ 



1 



i(zf), 



kzj l [1 + r z ] dz w 1 + r' z (zf). 



(35a) 
(35b) 



A deviation of r z from zero leads to a deviation of Zf from 
its unperturbed value. Using (|3"4"1) one can show that the 
shift is given by zj — Zf r z (zf), where Zf denotes 
value of the freeze-out temperature in the r z = case. 

Using equations (|3"4"|) and (j3"5")l we get for the asymp- 
totic value of the asymmetry 4 



T, 



~o l + r*(zf)-r' a (z f )+r' u (z f ) 

L i + (i + zf)n(zf) + faizf) 



(36) 



where the unperturbed value of the generated asymmetry 
in the strong washout regime is given by 



to _ 
1 l — 



2e N 2V27T 
7r 2 KZf exp(l) 



(37) 



Despite its simplicity, for n > 10 the analytical solution 
([56]) deviates from the exact numerical solution of 
by only a few percent and correctly describes the depen- 
dence of the asymmetry on r a and r z . For large k the 
inverse freeze-out temperature z f ^> 1 and the term pro- 
portional to Zf- 1 in the denominator can be neglected. 



The gauge invariance of r a and r± follows from the gauge invari- 
ancc of the (inverse) temperature and of <r and §| . 



4 On very small scales the thermodynamic quantities rapidly oscil- 
late and the solution obtained by the method of steepest descent 
under the assumption that the functions under the integral are 
smooth and do not oscillate is not valid. 
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B. The baryon— to— photon ratio 



transform of rv 5 



To complete the analysis we have to evaluate r a and 
Ti- Introducing the dimensionless temperature pertur- 
bation O = 4j£ one can relate z to its background value 
z = zq(1 + 6) _1 . Using the known dependence of 
zq on the time-coordinate in the radiation-dominated 
universe, zq = Hx° cx a, and the resulting relation 



zqZo 

dz 
IE 



H wc obtain for time derivative of 



H 



1 dz 



/g^dx z (l + 6)Vffoo 



which implies 



2Q + H~ 1 Q. 



(39) 



Expanding the exponent in (jlOp and neglecting second- 
order terms one can write a in the form 



a = a 3 (x°) (l- 3tp + J au], dx°) 



(40) 



where the scale factor a is related to the background tem- 
perature by a 3 cx T ~ 3 = T~ 3 (l + 36). Consequently 



r CT = 36 - 3tp + J au\dx° 



(41) 



Since r a and r-, are gauge invariant, it is sufficient 
to evaluate them in one particular gauge. In the lon- 
gitudinal gauge solution of the Einstein equations for 
the Fourier components of 6 and 6 in the radiation- 
dominated universe reads [l7| 



6 fc 



Cx 
2 

C 2 



SIM 



Co 



SIM 

X 

cos X 

X 



cos X 



SIM 



sin a; 



x 
sinx 



(42a) 



(42b) 



where x = ^= . The initial conditions are specified by C\ 
and C2, which are functions of k. The divergence of the 
macroscopic fluid velocity v = u l i can be calculated using 
the Einstein equation for the off-diagonal components of 
the Ricci tensor. In the absence of anisotropic stress, if) = 
<fr, and these simplify to (a</>) i = 4TrGa 2 (e + p )ui fvj\ . 
After use of the Friedmann equation, 3H 2 = 87rGa 2 £o, 
and some algebra we obtain for the corresponding Fourier 
components in the flat radiation-dominated universe 



2a- 



(acpk). 



(43) 



Trading the integration with respect to the time coor- 
dinate for the integration with respect to x and taking 
into account the ip = <f) equality we find for the Fourier 



r* k = 36 fc - 30 fc + / -^-—(x<t> k )dx = 0. (44) 



Using the relation 6 = i-5^6 = Hx^ we obtain a very 
simple expression for the Fourier transform of r-, 6 



b k + 26 fc + x- 



dx 



' 2 (Cx sin x + C2 cos x) 



(45) 

The lepton to photon number ratio, tjl = T^/T 7 , can 
now be easily read off from (|2"9"|) and (|3^|) . The quan- 
tity of interest is the departure of the asymmetry from 
its background value, r n = r\hj"rf L — 1. Since the baryon 
asymmetry is linearly related to the lepton one, r v also 
characterizes the baryon isocurvature perturbations. Ne- 
glecting the term proportional to zj 1 in the denominator 
of 



we find in the strong washout regime 7 
r n = r' u (z f ) - r*(z/) - 5^(2/). 



(46) 



The calculation of the CMB anisotropies is performed in 
momentum space, so that one is interested in the Fourier 
components of (|46j) . The Fourier transform must be eval- 
uated at the inverse freeze-out temperature Zf] it cor- 
responds to Xf = ZfX, where X = ^= is the value x 

takes at z = 1, and 1Z = Hr 1 cx MpiM^ 2 is the Hubble 
scale at this stage. The derivative of (|45|) with respect 
to the inverse temperature can be obtained upon use of 
the zero-order relation ^ = X-§^- Using the introduced 
notation we find for the Fourier transform of r v in the 
strong washout regime 

r Vk = X [£r u + l(2Cx z f - C 2 ) sxnXz f 
+ \{2C 2 z f + C 1 )cosXz f ] . 

In this regime the heavy neutrinos are tightly coupled to 
the SM species and it is natural to expect that r u van- 
ishes. The typical size of the perturbations is given by 
the Hubble radius 1Z. On larger scales r Vk is suppressed 



(47) 



5 Equation 1144 ( implies that for the photons T 7 = = is 
constant in the early universe, just as expected. Since photons 
are in equilibrium in the early universe, the right— hand side of the 
Boltzmann equation for photons vanishes and T 7 is conserved. 

6 A calculation using the solution of the Einstein equations for 
small perturbations in the synchronous gauge given in [23j leads 
to the same expressions for r a and r^; this confirms the gauge 
invariance of the introduced functions. 

7 Since the contribution of the heavy neutrinos to the energy den- 
sity is neglected, the baryon-to— photon ratio is related to the 
baryon-to-entropy ratio by a constant factor at this stage of 
the universe's evolution. Later, as the SM species decouple as 
they become massive, the baryon-to— entropy ratio remains con- 
stant, wheres the baryon-to— photon ratio changes. This change 
is absorbed into the definition of the background value of the 
asymmetry and does not modify the expressions for . 
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by the overall factor X, which becomes small. On su- 
perhorizon scales r Vk sa jX '{2C2Z f + C%) and vanishes as 

K — > OO. 

From (|29|) it follows that in the weak washout regime 



v dx 



exp(-ra 2 /2X 2 )dx , 



(48) 



where we have used the relation x = zX . That is, a 
nonzero r v is generated in this regime only if the velocity 
of the right-handed neutrino deviates from that of the 
SM species. The typical scale and amplitude of r u are 
model-dependent. Nevertheless one can make a plausible 
guess about these quantities: if perturbations in the SM 
and heavy neutrino components are created by the same 
mechanism, then the typical amplitude and scale of r u 
cannot exceed those of the velocity perturbation of the 
SM species. Since we are interested in large-scale per- 
turbations, which correspond to small X, it is sufficient 
to know the behaviour of r u in the vicinity of x — 



Q n X 



(49) 



If, for example, vanishes in the longitudinal coordi- 
nates 8 , then as follows from (|4*3")l o_i = 3C*2, ao = \C\, 
a,2 = 7C1 . . . , whereas if ity, vanishes in synchronous co- 
ordinates, then a\ = |C 2 , 04 = . . . (see Eq. (15) 
in (23|). To integrate the terms with negative n we cut 
the integral at some x m i n such that a„x™ in <C 1 to in- 
sure the applicability of the pcrturbative treatment. We 
furthermore consider only very large scales X <C x m in- 
Then 



r Vh = a^n^nX + 0(X 2 ) 



(50) 



That is, in the weak washout regime r u is also suppressed 
on large scales. 



C. BBN and CMB 

In the radiation-dominated universe the horizon size 
grows much faster than the scale factor, and conse- 
quently the generated perturbations are not protected 
from smoothing by baryon diffusion at later stages of the 
universe's evolution. 

Let us first consider the strong washout regime. Rela- 
tion (|47|) holds for arbitrary mass of the Major ana neu- 
trino. For the right-handed neutrino mass Mi ~ 10 9 
GeV, scales of the order of 1Z at the time of leptogenesis 
correspond to (roughly speaking) 10 -4 m at neutrino de- 
coupling. This is much smaller than the baryon diffusion 
length ~ 200 m [25| at this stage. Therefore in this case 



The analysis performed in 12411 in longitudinal gauge demon- 
strates that velocity perturbation of the heavy dark matter par- 
ticles are small compared to those of the massless species. 



the perturbations of the baryon to photon ratio are com- 
pletely smoothed out before nucleosynthesis starts. In 
the electroweak-scale resonant leptogenesis scenario [26| 
the Hubble scale at the time of leptogenesis corresponds 
to (roughly speaking) ~ 10 km at 0.1 MeV. At the latter 
temperature, which corresponds to the beginning of nu- 
cleosynthesis, the neutron diffusion length ~ 200 km. In 
other words, the baryon perturbations are not smoothed 
out on scales which correspond to X < 0.1 at the time 
of leptogenesis. Given the typical size of the integra- 
tion constants Ci^ ~ 10~ 5 — 10~ 4 the amplitude of the 
perturbations not washed out by baryon diffusion is too 
small to affect nucleosynthesis in an observable way. 

After the end of nucleosynthesis and before the begin- 
ning of recombination there are no neutral baryons and 
the diffusion length substantially decreases. Thus the 
baryon isocurvature perturbations which have not been 
smoothed out survive till recombination - the time of 
CMB formation. At this stage small scale perturbations 
blueof the photon gas undergo a mechanism called dif- 
fusion damping or Silk damping [27I . [28[ . The photons 
diffuse from regions with higher density to lower den- 
sity regions and thereby effectively smooth out baryon 
inhomogeneitics due to their large mean free path. This 
damps the amplitudes of acoustic oscillations by a factor 



D = exp 



Tdx { 



exp 



Me 

M 



2/3 



(51) 



where T is the damping rate, t° the present time and 
M c is called the critical mass that can only be obtained 
approximately from the integral over the damping rate. 
For a radiation-dominated universe during the phase of 
acoustic oscillations, one obtains 29] 



M c ks 8 x 10 13 (fi ^o)^ 1/2M (i 



(52) 



with density parameter f2o and rescaled Hubble param- 
eter ho evaluated at the present time and Mq denoting 
the solar mass. Sufficient damping is then obtained for 
(M c /M) of the order of 10 or larger. Translated into 
length scales, this means that all scales 



A D < 3 x 10 1 



(53) 



are affected by Silk damping at T ~ 10 eV, which roughly 
corresponds to the temperature at recombination, and 
only perturbations on scales larger than that can sur- 
vive. For the right-handed neutrino mass Mi ~ 10 9 
GeV, scales of the order of JZ at the time of leptogene- 
sis correspond to scales of the order of 10 m at recom- 
bination. This is obviously smaller than the damping 
scale Xd by so many orders of magnitude that no ob- 
servable imprint on the CMB can be produced. Turning 
the argumentation around, we can calculate the mass 
scale Mi that would be necessary to obtain perturba- 
tions on large enough scales to prevent them from being 
smoothed out by Silk damping, taking into account that 
1Z scales as T~ 2 in a radiation dominated universe. The 
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result is very disappointing: one would need a mass scale 
of Mi ^100 eV which is far below the scale where the 
sphalcrons arc in thermal equilibrium and can convert 
a lepton asymmetry into a baryon asymmetry. There- 
fore we conclude that in the strong-washout regime, the 
baryon inhomogeneities coming from the spatial fluctu- 
ation of the efficiency of leptogencsis can neither affect 
BBN nor the CMB in an observable way. 

Given the estimate (|50|) . the conclusions drawn in the 
case of the strong washout regime also apply to the case 
of the weak washout regime. 

IV. SUMMARY AND CONCLUSIONS 

In this short paper we have investigated the influence 
of the primordial perturbations of the energy density and 
metric on the efficiency of thermal leptogenesis. 

To perform the analysis we have derived an integral 
form of the Boltzmann equation for a gauge-invariant 
combination T of particle number density, the macro- 
scopic fluid velocity and space-time metric. In the ab- 
sence of the collision terms T is conserved to linear ap- 
proximation and replaces the particle number density per 
comoving volume in the inhomogeneous universe. 

Using the integral form of the Boltzmann equation ap- 
proximate analytical solutions for the lepton asymmetry 
in the weak and strong washout regimes have been found. 

In the strong washout regime the generated baryon 



isocurvature perturbations are correlated with the per- 
turbations of the energy density and metric. In the 
weak washout regime the solution deviates from the back- 
ground one only in the presence of the heavy neutrino 
velocity perturbation r u . The typical scale of the pertur- 
bations at the time of leptogenesis is set by the Hubble 
scale 7Z in both regimes. The amplitude is inversely pro- 
portional to the scale. 

The generated perturbations are smoothed by baryon 
diffusion at the later stages of the universe's evolution. 
The scale of the generated perturbations is comparable 
to the neutron diffusion length at the beginning of BBN 
only if the right-handed neutrino is as light as in the 
electroweak-scale resonant leptogenesis scenario. How- 
ever even in this extreme case the generated perturba- 
tions are smoothed out by the Silk damping at the re- 
combination. For this reason the influence of the baryon 
inhomogeneities coming from the spatial fluctuation of 
the efficiency of leptogenesis on BBN and CMB is unob- 
servably small. 
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